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LABORATORY  SIMULATIONS  OF  CONTROLLED  ENERGETIC 
ELECTRON-BEAM-PLASMA  INTERACTIONS  IN  SPACE 


L  INTRODUCTION 

From  the  time  of  the  early  experiments  involving  plasma 
1-4 

oscillations  there  has  been  continued  interest  in  the 
collective  interactions  of  charged-particle  beams  and  plasmas. 

This  interest  has  evolved  from  studies  of  various  small- 
scale  laboratory  gas  discharges  into  research  on  controlled 
thermonuclear  reactions,  and  recently  into  the  application 
of  charged-particle  beams  in  space.  These  efforts  have  made 
it  strikingly  clear  that  rarely  is  a  beam-plasma  system 
quiescent.  In  fact,  it  is  a  characteristic  of  beam-plasma 
systems  to  manifest  a  number  of  Instability  mechanisms  with 
attendant  and  varied  levels  of  plasma  turbulence. 

Within  the  space  science  community,  interest  in  artificially- 

injected  electron  beams  came  about  largely  from  a  perspective 

that  viewed  the  beams  as  diagnostic  devices  to  probe  various 

5-14 

aspects  of  ionospheric  and  magnetospher ic  phenomena 
The  applications  included  the  mapping  of  geomagnetic  field 
lines,  detection  of  geomagnetic  conjugates  by  the  generation 
of  artificial  aurora,  the  study  of  beam  spreading,  atmospheric 
excitation  and  ionization  processes  and  the  measurement  of 
magnetic  field-aligned  potentials.  (See  Ref.  14  for  an 
excellent  bibliography  and  chronology  of  rocket-borne  beam 
experiments  and  their  associated  results.)  On  the  other 
hand,  interest  in  nonlinear  beam-plasma  processes  arose  in  a 
less  direct  fashion.  The  early  rocket  experiments  resulted 
in  a  number  of  puzzles,  not  only  in  unexplained  beam-particle 
behavior  and  vehicle  neutralization  processes  but  also  in 
the  characteristics  of  emitted  plasma  waves.  In  some  cases 
the  results  were  at  such  variance  with  expectations  that  a 
beam-plasma  instability  process  was  suspected  of  playing  a 
dominant  role.  A  series  of  laboratory  simulations^-^  then 
accumulated  substantial  information  which  now  seems  to  have 
definitively  Identified  the  fundamental  process  as  the  beam- 
plasma-discharge^®-^  (BFD)  . 
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The  BPD  is  the  result  of  a  two-stream  beam-plasma 
instability  that  appears  at  a  critical  beam  current  I£ 
given  approximately  by 

I  [Amps]-6(10'11)V.  3/2/PLB0*7  (1) 

C  D 

where  the  beam  energy  is  in  volts,  the  superimposed 
magnetic  field  B  in  gauss,  the  pressure  P  in  torr,  and 
system  length  L  (gun-to-collector  distance)  in  meters* 

(Equation  (1)  is  an  empirical  fit  to  the  original  lab¬ 
oratory  studies17  for  parallel  injection  and  represents 
that  data  to  an  accuracy  better  than  40Z.)  This 
critical  current  level  yields  an  increased  ion-pair  production 
rate,  enhanced  3914&  emission,  a  modification  of  the  primary 
beam  velocity  and  emission  of  hf  waves. 

In  the  pre-BPD  case  (I  <  I  ) ,  the  beam  energy  distribution 

C  20  21 

is  sharply  peaked  at  the  acceleration  potential  *  ,  and 

the  hf  wave  emissions  are  usually  localized  to  narrow-bands 

slightly  greater  than  the  electron  cyclotron  frequency 

e  17 

o>c/2tt,  with  higher  harmonics  of  decreasing  amplitude 

The  amplitude  and  frequency  of  the  narrow-band  emission  is 

observed  to  increase  with  beam  current16.  Assuming  that 

density  varied  with  beam  current,  the  latter  result  suggested 

that  the  waves  were  at  the  upper  hybrid  frequency  oi  / 2 tt 

rather  than  at  / 2 tt  .  In  addition  to  the  narrow-band  emissions 

c 

there  are  low  frequency  waves  concentrated  at  the  low  end  of 

the  spectrum,  rising  steeply  below  Che  ion  plasma  frequency 

coi/2ir(^  100  kHz)  . 

P  > 

When  the  BPD  is  ignited  C I — I c ) »  the  entire  band  up  to 

the  electron  cyclotron  frequency  is  occupied  by  intense 

emissions  of  featureless  spectral  character is ics .  In  addition, 

an  intense  higher-frequency  band  of  emissions  appears  at  or 

near  the  local  plasma  frequency.  The  plasma  frequency  waves 

22 

are  most  intense  within  the  beam  column  (%  10  volts/meter), 
falling  off  sharply  from  the  beam  center  to  the  beam  edge 
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('V/  40  dB/meter)  .  Under  these  conditions,  the  beam  electrons 

are  redistributed  in  energy20’21  with  the  original  peak  at 

the  accelerating  potential  (typically  800  volts  to  2  keV) , 

broadened  10-15%  FWHM  (full-width-half-max)  with  an  enhanced 

low  energy  tail  in  the  E<  100  eV  domain. 

Following  the  characterization  of  the  BPD  by  wave 

emissions,  beam  energy  redistribution  and  the  various  parametric 

dependencies,  space-simulation  efforts12  then  focussed  on 

determining  a  plasma-density-dependent  criterion  for  BPD 

17  18 

ignition.  Early  thoughts  ’  suggested  that  the  BPD  was 

0 

ignited  when  the  local  electron  plasma  frequency  was 
greater  than  or  equal  to  the  local  electron  cyclotron  frequency 

<oe  (i.e.,  coe-  coe)  .  This  condition  was  found  to  be  consistent 

c  P  c  23  24 

with  the  most  recent  results  ’  which  established 


w  -  5.4  oie  ... 

P  p  (2) 

as  the  ignition  criterion  appropriate  to  the  laboratory 
efforts  at  the  Johnson  Space  Center.  The  density-dependent 
criterion  also  agreed  with  theoretical  calculations  which 
considered  the  BPD  to  be  produced  by  large  amplitude  electron 
plasma  waves  resulting  from  the  beam-plasma  interaction. 

While  many  of  the  laboratory  results  were  finding 

23  25 

acceptable  theoretical  descriptions  ’  ,  additional  work 

was  required  to  more  completely  define  the  basic  plasma 
processes  and  plasma  electron  distributions.  Critical 
issues  included  direct  measurements  of  enhanced  density 
levels  and,  more  Importantly,  a  definitive  experimental 
determination  of  suprathermal  electron  populations  at  BPD 
threshold.  Until  recently  the  latter  issue  had  been  mostly 
a  matter  of  theoretical  speculation.  With  this  in  mind  and 


with  interest  in  exploring  spectral  energy  transfer  processes 
and  the  conversion  of  one  type  of  fluctuation  into  another, 
we  have  conducted  a  number  of  follow-on  experiments.  These 
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experiments  dealt  primarily  with  the  plasma  electron  populatio 
its  density  and  energy  distributions,  and  levels  of  turbulence 
with  respect  to  the  various  BPD  states  (pre-,  threshold-, 
and  solid  BPD).  The  experimental  configuration,  results  and 
extrapolations  to  space-borne  applications  are  treated  in 
subsequent  sections. 

IL  EXPERIMENT  CONFIGURATION  AND  DIAGNOSTIC  TECHNIQUE 

The  experiment  was  conducted  in  the  large  vacuum  chamber 

facility  at  the  NASA  Johnson  Space  Center  in  Houston,  Texas 

with  a  configuration  (Figure  1)  similar  to  that  employed  in 

the  earlier  space  simulations^  .  The  chamber,  with  base 

pressures  'V*  (10~^)  torr,  was  equipped  with  large  current- 

carrying  coils  to  generate  magnetic  fields  up  to  2.1  gauss. 

A  steerable  Pierce-type  electron  gun  with  a  perveance  of 
—  6  —3/2 

^  1.4(10  )AV^  was  mounted  near  the  chamber  floor  on  a 

movable  cart  that  allowed  the  beam  to  be  injected  upward  and 
parallel  to  the  magnetic  field  B,  and  to  be  terminated  on  a 
gridded  3  x  3m  collector  suspended  20  m  above  the  gun 
aperture . 

The  chamber  was  also  equipped  with  a  position-controlled 

29  3 

cylindrical  pulsed-plasma-probe  system  (acronym  P  )  that 
could  be  continuously  varied  in  its  radial  position  relative 
to  the  beam  center.  All  radial  traversals  were  along  the 
local  magnetic  meridian  at  a  height  of  8  m  above  the  gun 
aperture.  Care  was  taken  to  guarantee  that  radial  profile 
information  was  not  distorted  by  magnetic-aspect  sensitivities 

3 

The  P  technique  provided  measurements  of  relative 

electron  density  on  a  1  msec  time  scale  while  simultaneously 

generating  a  "conventional"  Langmuir  probe  characteristic 

for  determination  of  absolute  N  ,  T  and  plasma  potential 
3  e  e 

V  .  In  the  P  procedure  a  chain  of  voltage  pulses  are  applied 


to  the  probe*  The  pulse  amplitudes  follow  a  sawtooth 
envelope  and  generate  the  data  pairs  (I,V)  for  the  conventional 
Langmuir  probe  I-V  characteristic.  During  the  interpulse 
period  the  probe  is  held  at  a  fixed-baseline  level,  Vg ,  in 
the  electron-saturation  portion  of  the  characteristic.  The 
running  measurements  of  Ig,  during  the  baseline  period  Vg, 
provide  a  measure  of  relative  N  variations  (assuming  1  - 

.  6  O 

g  a  ^ 

I  «  N and  the  fundamental  data  set  for  further  analysis 
of  density  fluctuation  power  spectra.  The  procedure  also 
mitigates  hysteresis  effects  known  to  degrade  conventional 
Langmuir  probe  procedures.  (References  24  and  29  and 
associated  bibliographies  are  recommended  for  more  detailed 

3 

information  on  the  P  technique.)  The  probe  was  fabricated 
from  0.8  mm  diameter  tungsten  wire,  20  cm  in  length.  It  was 
separated  from  an  extension  arm  by  an  2.5  cm  long  coaxial 
guard  electrode  that  was  driven  at  the  same  potentials  as 
the  probe.  The  plasma  currents  collected  by  the  probe  were 
measured  by  an  electrometer  system  that  provided  linear  bi¬ 
polar  outputs  over  a  current  range  covering  five  orders  of 
magnitude . 

m.  EXPERIMENTAL  RESULTS 
Bulk  Plasma  Features 

To  establish  a  framework  for  the  plasma  properties  to 

be  discussed  here  and  in  subsequent  sections  it  is  useful  to 

adopt  the  three  stage  description  of  BPD-buildup  proposed  in 

1 8 

the  work  of  Getty  and  Smullin  . 

Stage  1  (Pre-BPD) :  The  beam  behaves  in  the  classical 
sense  of  single-particle  interactions,  creating  ion-electron 
pairs  upon  collision  with  ambient  neutrals.  The  beam- 
created  plasma  is  low  in  density  and  increases  linearly  with 
beam  current. 

Stage  2  (Approaching  Threshold) ;  Collective  effects 
begin  to  take  place  with  the  beam-plasma  system  generating 
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oscillations  that  first  appear  near  the  electron  cyclotron 
frequency.  These  oscillations  (may)  produce  plasma  elec¬ 
trons  in  excess  of  that  produced  directly  by  the  beam. 

Stage  3  (BPD) :  More  intense  oscillations  are  generated 
at  and  above  the  electron  plasma  frequency.  This  is  believed 
to  be  the  dominant  mode  of  interaction,  resulting  n  substantial 
increases  in  plasma  production. 

While  this  staging  process  was  developed  to  describe 

the  temporal  buildup  for  a  pulsed  beam  operation  with  I  > 

Ic ,  it  has  merit  for  the  description  of  dc  discharges  and 

the  development  toward  BPD  at  various  beam  current  levels. 

In  Figure  2  are  plotted  maximum  values  of  plasma  density 

measured  along  radial  plasma  profiles  for  fixed  beam  energy 

(1.3  keV)  and  selected  levels  of  beam  current.  The  data  are 

seen  to  support  a  threshold  concept  in  which  the  rate  of 

plasma  production  is  enhanced  at  some  critical  current  level 

(12  ma  in  Fig.  2).  While  the  enhanced  ion-pair  production 

rate  may  be  characterized  by  an  avalanche  break-down  as 

suggested  in  some  of  the  earlier  space-simulation  efforts, 

there  are  cases  in  which  the  transition  has  been  found  to  be 

more  gradual.  Experience  accumulated  to  date  indicates  that 

at  threshold  N  (I^I  ) /N  (I^I  )  never  exceeds  20,  with  factors 
e  c  e  c 

in  the  range  of  2-5  being  more  typical. 

Radial  profiles  of  plasma  density  for  pre-,  threshold-, 
and  solid-BPD  conditions  are  presented  in  Figures  3  and  4 
for  ■  1.9  and  1.3  keV,  respectively.  The  abscissa  is 
time  relative  to  the  start  of  each  radial  traversal  and  the 
o  inate  is  the  logarithm  of  relative  electron  density  as 
determined  by  baseline  electron-saturation  currents 
collected  by  the  cylindrical  pul sed-plasma-probe .  At  the 
start  of  each  traversal  the  probe  was  at  its  outermost 
position  relative  to  the  center  of  the  chamber  and  the  beam 
axis.  As  time  increased  the  probe  was  moved  into  and 


through  the  beam  at  a  rate  ^7.5  cm/sec;  after  penetration 
through  the  beam  center  the  traversal  system  was  reversed, 
allowing  a  second  measurement  of  the  density  profile  as  the 
probe  moved  back  to  its  original  outermost  position.  With 
this  procedure  the  probe's  maximum  penetration  into  and 
through  the  beam  is  identified  by  the  symmetry  point  (see 
e.g.,  top  panel  of  Figure  3)  in  the  "double"  profile. 

The  top  panel  in  Figure  3  is  typical  of  a  pre-BPD 
profile  in  that  the  beam  center  is  well  defined  and  of 
relatively  narrow  cross  section.  At  BPD  threshold  the 
plasma  expands  in  cross  section  generally  eliminating  the 
narrow  beam  profile  characteristic  of  the  pre-BPD  state.  As 
the  beam  current  is  increased  above  the  threshold  level,  the 
plasma  density  increases  and  the  radial  profiles  take  on 
various  configurations.  For  example,  in  the  bottom  panel  of 
Figure  3,  the  center  is  depleted  with  an  increased  density 
level  at  its  edges.  This  configuration  has  been  observed  in 
a  number  of  cases  and  is  thought  to  be  the  result  of  electron 
heating  in  the  beam  core  with  generally  higher  diffusion 
rates  and  associated  losses.  The  fact  that  a  similar  depletion 
is  not  observed  in  the  bottom  panel  of  Figure  4  indicates 
that  our  understanding  is  less  than  complete.  In  fact,  the 
results  in  Figure  4  have  been  selected  to  present  information 
at  a  higher  magnetic  field  (1.5  gauss  in  Fig.  4  and  0.9 
gauss  in  Fig.  3)  and  to  illustrate  that  the  description  of 
beam  plasma  profiles  offered  in  connection  with  Figure  3 
should  be  considered  nominally  in  agreement  with  overall 
observations  but  not  universally  applicable  to  the  various 
levels  of  the  BPD.  (For  purposes  of  perspective  in  the 
development  of  BPD,  it  is  noted  that  the  profiles  in  Figures 
4,  upper,  middle  and  bottom  panels,  correspond  to  the  8,  12 
and  38  ma  cases,  respectively  in  Figure  2.) 
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De natty  Fhictuatioiu 

The  density  profiles  shown  in  Figures  3  and  4  are 
constructed  of  one-second  averages  along  the  probe  traversals. 
The  actual  profiles  are  not  as  uniform  as  the  figures  would 
suggest,  but  are  subject  to  rather  large-scale  radially- 
dependent  density  fluctuations.  As  an  illustration.  Figure 
5  presents  half-second  segments  of  density  fluctuations  at 
four  positions  relative  to  the  beam  center  under  the  pre-BPD 
conditions  corresponding  to  the  top  panel  in  Figure  3.  In 
addition  to  frequency  components  easily  discernable  on  a 
half-second  time  frame,  the  data  in  Figure  3  show  that  the 
relative  fluctuation  amplitude  increases  with  radial  distance 
from  the  beam  center.  The  frequency  components  have  been 
fast-Fourier  analyzed  (see  Fig.  6)  to  determine  the  spectral 
distribution.  The  spectrum  shows  a  dominance  at  the  low 

frequency  end  (1  Hz  ^  f  ^  20  Hz),  with  components  at  higher 

'  -3 

frequencies  (f  -►  250  Hz)  exhibiting  an  f  dependence. 

The  radial-dependence  of  fluctuation  amplitudes  is 
presented  from  two  perspectives  in  Figures  7  and  8.  Both 
figures  cover  the  pre-BPD  conditions  in  the  top  panel  of 
Figure  3.  In  Figure  7  is  plotted  the  maximum  spectral 
intensity  measured  in  a  radially-dependent  series  of  fast- 
Fourier-analyzed  density  fluctuations.  Each  data  point 
represents  a  Fourier  analyzed  block  of  density  irregularities 
at  specific  positions  in  the  profile  and  the  variations  from 
point-to-point  reflect  the  nondeterministic  nature  of  the 
plasma  turbulence.  The  solid  curve  is  a  smooth  fit  to  the 
results  and  defines  more  quantitatively  the  radially-dependent 
fluctuation-amplitude  distribution  evident  in  the  raw  data 
(Figure  5) .  The  results  show  that  the  relative  spectral 
intensity  of  density  fluctuations  in  the  beam-plasma  system 
is  a  factor  of  10-20  times  greater  at  a  3  meter  distance 
from  the  beam  center  than  it  is  within  the  beam  core  itself. 
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A  similar  observation  (though  somewhat  different  on  a 
quantitative  scale)  is  reflected  in  the  results  of  Figure  8 
where  |<5N  I  /N  («|5l_|/I_)  is  plotted  in  a  horizontal  frame 
that  corresponds  identically  with  the  abscissa  (and  the 
conditions)  in  the  top  panel  of  Figure  3.  The  plotted 
points,  which  are  absolute  deviations  from  sliding  linear 
detrends,  show  that  the  fluctuations  are  largest  outside  the 
beam  center  with  the  greatest  frequency  of  occurrence  at 
fluctuation  amplitudes  -  1.0  (i.e.,  less  than  or  equal  to 
100%  fluctuations).  On  a  less  frequent  basis  the  amplitude 
fluctuations  can  extend  to  2.0  (i.e.,  200%).  Toward  the 
beam  center  the  fluctuation  amplitude  is  observed  at  significantly 
reduced  levels,  corresponding  approximately  to  0.2  (20%). 

The  trend  in  density  fluctuations  represented  in 
Figures  5,  7  and  8  can  be  considered  typical  for  most  beam- 
plasma  interactions  studied  to  date  in  the  Johnson  Space 
Center  facility. 

In  an  effort  to  identify  the  wave  modes  responsible  for 
the  observed  density  fluctuations,  it  is  of  benefit  to 
return  to  the  subject  of  spectral  distributions.  In  Figure  9 
is  presented  a  fast-Fourier  analyzed  sample  of  density 
fluctuations  observed  under  the  threshold  conditions  in  the 
middle  panel  of  Figure  4.  This  sample  has  been  selected 
because  it  reflects  substantial  contributions  from  higher 
frequency  components  (f  @  55  Hz)  and  because  additional 
attention  will  be  given  to  that  same  threshold  condition  in 
the  next  section.  In  composite,  the  results  in  Figures  6 
and  9  are  typical  of  the  spectral  observations  accumulated 
to  date:  (1)  Generally  the  density  fluctuations  are  dominated 
by  low  frequency  components  ('V'  1-50  Hz)  with  relatively 
constant  levels  of  Intensity;  (2)  While  observed  frequencies 
extend  to  250  Hz  (higher  frequency  components  have  as  yet 
not  been  studied) ,  the  region  at  or  near  the  ion  cyclotron 
frequency  (50  Hz  in  Figure  9)  appears  to  manifest  the 
highest  frequency  components  of  substantial  spectral  intensity; 
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(3)  When  observed,  the  ion  cyclotron  coaponent  is  never  a 
narrow  peak  in  the  spectral  distribution  but  rather  a  broad 
feature  centering  at  or  near  <u*/2ir. 

While  the  current  set  of  observations  does  not  provide 
the  wave  number  measurements  to  complement  the  frequency 
spectrum  results,  the  domain  of  contributing  wave  phenomena 
can  be  isolated  by  reference  to  Figure  10.  For  parametric 
ranges  encompassing  the  conditions  in  the  space-simulation 
experiments,  the  figure  presents  the  results  of  linear 
theory  for  the  wavenumber  dependence  of  low  frequency 
fluctuations  extending  from  the  lowest  frequency  ion  acoustic 
waves  (at  1  Hz)  to  whistler  and  Alfven  waves  at  100  kHz. 

Limiting  the  discussion  to  the  more  intense  frequencies 
in  the  measurements  of  density  fluctuations  (i.e.,  1  Hz  - 
f  ^  55  Hz)  and  to  wavelengths  compatible  with  the  system 
size  (20  meters  and  less) ,  it  appears  that  there  are  two 
likely  contributions: 

(i)  the  longitudinal  ion  cyclotron  waves  branching 
into  the  ion  acoustic  mode,  and 

(ii)  the  lower  branch  of  the  drift  wave  mode  which  in 
itself  covers  all  frequencies  in  the  current  set  of  observa¬ 
tions. 

The  measurements  support  both  possibilities,  but  the 
existence  of  gradients  in  density  and  temperature  (detailed 
below)  favor  the  drift  wave  mode. 

A  Complete  Plasma  Profile 

Thus  far,  attention  has  been  directed  to  absolute 
densities,  density  fluctuation  power  spectra  and  associated 
macroscopic  features  relevant  to  the  various  stages  of  beam- 
plasma  interactions.  However,  a  more  complete  analysis  of 
the  interactive  beam-plasma  processes  requires  a  study  of 
the  possible  effects  on  electron  energy  distribution  functions 
Since  work  on  the  beam  electron  distribution  has  already 
been  reviewed,  we  concentrate  here  on  the  thermal  electron 
population  and  the  existence  of  suprathermal  electron  tails 


I 


as  suggested  in  recent  theoretical  descriptions  of  the  beam- 
23 

plasma-discharge  .  A  composite  experimental  profile  of 

pertinent  plasma  parameters  measured  simultaneously  by  the 
3 

F  diagnostic  system  is  presented  in  Figure  11  for  the 
threshold  condition  corresponding  to  the  middle  panel  of 
Figure  4.  All  data  for  the  composite  profile  were  extracted 

3 

from  P  characteristics,  with  a  sample  in  Figure  12  showing 
a  two-component  electron  population  within  the  beam  core. 

All  related  analyses  assume  both  components  to  be  Maxwellian 
and  characterized  by  an  electron  temperature.  Accordingly, 
the  retarding-f ield  and  plasma  potential  domains  were  analyzed 
with  conventional  procedures  normally  employed  in  such 
environments . 

The  results  show  the  suprathermal  electron  population 

confined  to  a  region  within  1  meter  of  the  beam  center,  with 

the  total  population  falling  to  its  e-fold  level  in  60-70 

cm.  The  steep  radial  dependence  of  the  suprathermal  component 

is  similarly  reflected  in  the  relative  profile  of  the  total 

5  -3 

electron  density.  At  its  peak  density  level,  8  (10  )cm  , 

the  suprathermal  component  is  20%  of  the  total  population. 

The  temperature  of  the  suprathermal  component  is  seen 

to  rise  very  sharply  at  the  first  appearance  of  the  hot 

electrons  near  1.3  meters  and  then  levels  off  to  a  near 
4  o 

constant  5.3  (10  )  K  level  near  the  center.  The  cold  electron 

A  o 

component  is  1.5  (10  )  K  at  its  peak,  falling  off 

exponetially  with  radius,  in  agreement  with  previously 

2  6 

reported  BPD  observations 

The  figure  also  reveals  an  interesting  feature  relative 
to  the  coincident  measurements  of  density  fluctuations.  The 
dashed  line  marking  the  spatial  onset  of  density  fluctuations 
(an  exact  1:1  temporal  correspondence  can  be  located  in  the 
time-dependent  radial  traversal  in  the  central  panel  of 
Figure  4)  is  co-located  with  an  exponential-like  increase  in 
the  relative  spectral  intensity  and,  more  significantly, 
with  the  steep  temperature  gradient  in  the  suprathermal 
electron  population.  This  observation  suggests  that  the 
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turbulence  is  initiated  at  this  suprathermal  temperature 
boundary  and  expands  radially  with  increasing  amplitude  as 
the  waves  penetrate  lower  plasma  density  domains. 

One  final  observation  is  included  in  the  upper  panel  of 

Figure  11,  where  measurements  of  plasma  potential  are 

plotted  in  the  same  radial  format  as  density,  temperature 

and  spectral  intensity.  It  is  seen  that  the  plasma  potential 

is  a  well  behaved  function,  increasing  in  magnitude  (relative 

to  the  chamber  wall)  with  decreasing  radius.  The  overall 

potential  drop  is  only  0.6  volts,  corresponding  to  'V*  1  kT 

3  o  * 

for  an  average  cold  electron  temperature  of  8(10  ;°K. 

IV.  COMMENTS  AND  CONCLUSIONS 

Although  the  beam-plasma-discharge  has  had  a  relatively 
long  history  of  laboratory  investigations,  recent  space- 
related  efforts  have  brought  about  a  renewed  interest  with 
concerted  theoretical  and  experimental  work  focussing  on  the 
problem.  Accumulated  information  has  helped  expand  many  of 
the  original  concepts  and  several  new  quantitative  aspects 
of  the  problem  have  evolved.  The  progress  has  resulted  in 
part  from  the  improved  understanding  of  the  energization 
processes  of  the  ambient  plasma  in  its  Interaction  with  the 
energetic  electron  beam.  Another  major  contribution  has 
come  from  the  team  effort  in  attacking  the  problem  with  a 
perspective  that  concentrates  on  a  more  comprehensive  description 
(both  theoretical  and  experimental)  of  the  Interaction 
process,  including  diagnostics  of  the  waves,  fields,  particle 
densities  and  energy  distribution  functions.  While  the 
efforts  reported  in  the  present  investigation  have  dealt 
only  with  the  plasma-electron  properties,  an  understanding 
of  the  observations  was  made  possible  only  within  the  context 
of  all  accumulated  information. 

The  new  results  include: 

( i)  Quantification  of  density  enhancement  in  the  transition 
from  single-particle  to  the  collective  BPD  process.  While 
the  BPD-enhanced  ion-pair-production  rate  may  be  characterized 
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by  an  avalanche  breakdown  as  suggested  in  the  results  of 
Figure  2,  there  are  cases  in  which  the  transition  is  sore 
gradual.  Experience  accumulated  to  date  Indicates  that  at 
BPO-threshold  the  enhancement  never  exceeds  20,  with  factors 
in  the  range  2-5  appearing  more  typical. 

( ii)  Radial  profile  characterization  as  a  function  of 
pre-BPD,  threshold,  and  BPD  conditions.  The  pre-BPD  plasma 
electron  profile  is  well  defined  and  of  relatively  narrow 
cross  section.  The  beam-width  itself  is  of  the  order  tens 
of  centimeters  with  the  plasma  density  falling  exponentially 
to  the  chamber  walls.  At  BPD  threshold  the  plasma  column 
expands  in  cross  section  (1-2  meters),  generally  eliminating 
the  narrow  beam  profile  characteristic  of  the  pre-BPD  state. 

As  the  beam  current  is  increased  above  BPD-threshold ,  the 
plasma  density  increases  and  the  density  profiles  take  on 
various  configurations  with  evidence  for  electron  heating  in 
the  beam  core  and  generally  higher  diffusion  rates.  (See 
Figs.  3  and  4.) 

(iii)  Amplitude  and  spectral  classification  of  density 
fluctuations .  Density  fluctuations  are  dominated  by  low 
frequency  components  ('v*  1-50  Hz)  with  amplitudes  increasing 
with  distance  from  the  beam  center.  Fluctuation  levels  can 
reach  100%.  Candidate  wave  modes  Include:  (a)  longitudinal 
ion  cyclotron  waves  branching  into  ion  acoustic  modes,  and 
(b)  the  lower  branch  of  the  drift  wave  mode.  The  presence 
of  steep  gradients  in  density  and  temperature  favor  the 
latter  mode  as  dominant. 

( iv)  Identification  and  mapping  of  a  suprathermal  electron 

component  at  and  above  BPD  threshold.  Copious  amounts  of 

suprathermal  electrons  were  found  to  coincide  with  the  onset 

5  -3 

of  BPD.  These  electrons  (Ng  -  8(10  ) cm  at  5  eV)  are  part 

of  a  broad  distribution  which  extends  up  through  the  50-100  eV 
31 

range  .  They  result  from  strong  wave-particle  interactions 
and  trigger  the  enhanced  ionization  process  characteristic 
of  the  BPD.  Their  observation  represents  an  important  step 
in  experimental  studies  of  beam-plasma  interactions  and 
significantly  improve  the  developing  model  of  the  BPD. 
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ELECTRON  DENSITY  [107  cm  3| 


Fig.  2  —  Maximum  values  of  electron  density 
observed  in  radial  plasma  profiles  at  levels  of  beam 
current  spanning  the  pre-,  threshold-  and  solid 
BPD  regimes. 
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Fig.  4  —  Radial  profiles  of  plasma  density  for 
pre-  (upper  panel),  threshold  (middle  panel)  and 
solid  (bottom  panel)  BPD  conditions  at  Vb  ■ 
1.3  keV  and  B  *  1.5  gauss. 
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RUN  NUMBER  46. 
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Fig.  5  —  A  series  of  0.5  second  data  segments 
showing  density  fluctuations  at  four  positions 
within  the  radial  profile  shown  in  the  upper  panel 
of  Figure  3.  The  times  displayed  within  each  panel 
can  be  directly  cross  referenced  with  the  traversal 
time  (and  relative  position  with  respect  to  the  beam 
center)  in  the  top  profile  in  Figure  3. 
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POWER  SPECTRUM 


Fig.  6  —  FFT  analysis  of  a  2-second  segment  in  the 
profile  shown  in  the  upper  panel  of  Figure  3. 
The  Fourier  analyzed  data  sample  includes  the 
20.0-20.5  sec  segment  displayed  in  Figure  5. 
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RADIAL  DISTANCE  [Meters] 

Fig.  7  —  Radial  dependence  of  spectral  intensity  associated  with  the 
density  fluctuations  measured  under  the  conditions  of  Figure  3,  top 
panel  (pre-BPD  condition).  In  each  case  the  intensities  represent 
peak  values  in  the  overall  spectral  distribution. 
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DENSITY  FLUCTUATIONS 
|dN.|/N, 


TRAVERSAL  TIME  [sec]—** 

Fig.  8  —  Radial  distribution  of  density  fluctuations  l5Ne  l/Ne  occurring  in  the 
pre-BPD  conditions  of  Figure  3,  top  panel.  The  light  and  dark  vertical  stratifi¬ 
cations  are  a  result  of  computer-plotter  anomalies. 


FREQUENCY  (Hz) 

Fig.  9  —  FFT  analysis  of  a  2-second  segment  in  the  threshold-BPD 
profile  shown  in  the  middle  panel  of  Figure  4 
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Fig.  10  —  Dispersion  relations  for  low  frequency  fluctuations  in  a 
magnetoplasma.  The  following  parametric  values  have  been  as¬ 
sumed:  ne  -  nj  ■  106  cm-3 ,  T,  *  3000°K,  Mj  *  30  amu  and  B  - 1 
gauss. 
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Fig.  11  —  Composite  profile  of  measured  plasma  parameters  under  the  BPD-threshold 
conditions  in  the  middle  panel  of  Figure  4.  Actual  data  points  are  not  included  in  the 
figure  but  selected  error  bars  have  been  utilized  to  represent  the  range  of  inaccuracies 
in  the  smoothed  fit.  Note  that  all  curves  are  plotted  on  a  linear  scale. 
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Fig.  12  —  Pulsed -plasma-probe  characteristic  showing 
evidence  of  a  suprathermal  electron  component 
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